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^ ■ Abstract 

An extended version of the antisymmetrized molecular dynamics to study structure of p-sd shell 
hypernuclei is developed. By using an effective AN interaction, we investigate energy curves of ^Be, 
^ 3 C and ?°' 21 Ne as function of nuclear quadrupole deformation. Change of nuclear deformation 

a: 

caused by A particle is discussed. It is found that the A in p-wave enhances nuclear deformation, 
while that in s-wave reduces it. This effect is most prominent in ^ 3 C. The possibility of the parity 
inversion in ^°Ne is also examined. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



One of the unique and interesting aspects of hypernuclei is the structure change caused 
by the hyperon as an impurity. Many theoretical works have suggested such phenomena 
caused by A particle as the change of deformation [l-5], the shrinkage of the inter-cluster 
distance l|, |2j and the super-symmetric (genuine) hypernuclear state |2J, |6j-l9] . Owing to the 
experimental developments some of them have been observed in light p-shell hypernuclei. 



As examples, we can refer to the reduction of B(E2) in .1 
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r] and the identification of the 



super-symmetric (genuine) hypernuclear state in \Be 

Today, we can expect that a new experimental facility of Japan Proton Accelerator Re- 
search Complex (J- PARC) will reveal the spectral information of p-sd shell and neutron-rich 
hypernuclei. Since these normal nuclei have a variety of structure such as coexistence of shell 
and cluster structure 
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16( | and novel exotic clustering [17H21I]. there must be many in- 



teresting phenomena peculiar to hypernuclei to be found. Indeed, several pioneering works 
predicted exotic structure in sd-shell hypernuclei such as the parity inversion in ^°Ne [?] and 



various rotational bands in ^Ne 



221 ] . These works are based on rather limited knowledge 



on the AN interaction. Since our knowledge of AN interaction has been greatly increased by 

I3I . 23 -27 1, we are now able to perform more 



the recent theoretical and experimental efforts 13, 23 
quantitative and systematic study of the structure change in A hypernuclei. It will reveal 
the dynamics and many interesting aspects of baryon many-body problem. 

To perform systematic and quantitative study of s<i-shell and neutron-rich A hypernu- 
clei, we develop an extended version of the antisymmetrized molecular dynamics (AMD) 
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28H31j. which we shall call HyperAMD. AMD has been successful to describe various 



exotic structure of neutron-rich nuclei and highly excited states of stable nuclei. Therefore, 
HyperAMD is suitable to describe the structure change and exotic structure of hypernuclei. 

In this study, we introduce HyperAMD and focus on the change of nuclear quadrupole 
deformation caused by a A particle. By applying HyperAMD to ^Be, \ 3 C, ^°Ne and ^Ne 
with YNG-ND AN interaction [32] ]. it is found that A particle changes nuclear quadrupole 
deformation. While the A particle in s-wave reduces quadrupole deformation as expected, 
that in p-orbital increases it. Among the calculated hypernuclei, jfC has shown the most 
drastic change of the nuclear deformation. It is also found that the binding energy of A 
particle depends on the structure of the core nucleus. Namely, the A in s-wave coupled 
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to the deformed core nucleus has smaller binding than that coupled to spherical core. On 
the contrary, the A in p-wave coupled to the deformed core has larger binding than that 
coupled to the spherical core. This contradicts to the preceding study [7] in which the parity 
inversion of ^°Ne is predicted. 

This paper is organized as follows. In the next section, we explain the theoretical frame- 
work of HyperAMD. In the section III, the change of energy curves as function of quadrupole 
deformation are presented. The trend of the change and its origin are discussed. The final 
section summarizes this work. 



II. THEORETICAL FRAMEWORK 



In this section, we introduce the theoretical framework of HyperAMD. Compared to the 



coupled channel AMD which describes the multi-strangeness system [33], it has better de- 
scription of the hyperon single particle wave function, but it does not treat multi-strangeness 
and is limited to single A hypernuclei. 

A. Wave function 

A single A hypernucleus consists of A nucleons and a A particle is described by the wave 
function that is the eigenstate of the parity, 

*± = P^ int , (1) 

where P ± is the parity projector, and the intrinsic wave function \l/ int is given by the direct 
product of the A single particle wave function and the wave function of A nucleons \&at, 

*int = <PA®^N- (2) 

The nuclear part is described by a Slater determinant of nucleon single particle wave packets, 

* N = -^dettyifo)}, (3) 
M r j) = <M r j) -Xi-Vh ( 4 ) 

a=x,y,z ^ ' *■ * 

Xi = aiX^ + PiXh (6) 



r\i = proton or neutron, (7) 



where ipi is the i-th nucleoli single-particle wave packet consists of spatial spin \i an d 
isospin rji parts. The spatial part 4>i is represented by a deformed Gaussian. Its centroid Zi 
is a complex valued three-dimensional vector. The width parameters v a are real numbers 
and take independent value for each direction, but are common to all nucleons. The spin 
part is parameterized by the complex parameters aij and fa, and the isospin part is fixed to 
proton or neutron. Zi, v a , ccj and fa are the variational parameters of the nuclear part. 

To describe various wave functions of the A particle, the A single particle wave function 
is represented by the superposition of Gaussian wave packets, 

M 

V'aO) = ^2c m ip m (r), <p m (r) = 4>m{r) ■ Xm, (8) 



m=l 



(t>m{r) = Yl (~~~) expl-u a p(r - z m )l L (9) 

a=x,y,z \ / v 

Xm = a m Xt + b m Xh ( 10 ) 

p - ^. (11) 

m N 

Again, each wave packet is parametrized by the centroid of Gaussian z m , the spin direction 
a m and b m . z m , a m , b m and c m are the variational parameters of the A part. The width 
parameter v a are equal to those of nuclear part. The number of the superposition M is 
taken to be large enough to have the energy convergence of the variational calculation. 



B. Hamiltonian 



The Hamiltonian used in this study is given as 

H = f N + V NN + Vcoui + T A + Van - f g . (12) 

Here, T N , Ta and T g are the kinetic energies of nucleons, A particle and the center-of-mass 
motion. Since we have superposed Gaussian wave packets to describe the A single particle 
wave function, it is rather time consuming to remove the spurious motion of the center-of- 
mass exactly. To reduce the spurious motion, we keep the center-of-mass of wave packets at 
the origin of the coordinate, 

A M 

J2z t + J2vp~ z ™ = 0- (is) 

i=l m=l 
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We expect that the spurious energy is not large in s<i-shell hypernuclei, since the number 
of nucleons is much larger than the hyperon. A similar method is also applied in the other 
AMD studies [28 1 . 



Our model wave function is designed to describe the low-momentum phenomena as in the 
case of the conventional shell model and we shall use the low-momentum effective interaction. 



We have used the Gogny D1S interaction [34j as an effective nucleon-nucleon interaction Vnn, 
that has been successfully applied to the stable and unstable normal nuclei. As an effective 
AN interaction, we have used the central part of the YNG-ND interaction [32] . The YNG-ND 
interaction depends on the nuclear Fermi momentum kp through the density-dependence of 
the G-matrix in nuclear medium. In this work, we apply respectively fcp=1.14 and 1.17 fm _1 
for ^Be and \ 3 C, that are so determined to approximately reproduce the binding energy of A 
in s-wave. For ^°Ne and /^Ne, we apply the same value as ^ 3 C, since there is no experimental 
data. The Coulomb interaction is approximated by the sum of seven Gaussians. 

C. Frictional cooling method with constraints 

Using the frictional cooling method, the variational parameters of the model wave func- 
tion are so determined that the total energy is minimized under the constraints. We have 
imposed two constraints on the variational calculation. The first is on the nuclear quadrupole 
deformation parameter (3 that is achieved by adding the parabolic potential, 

(V )= v ^-fo)\ (14) 

to the total energy. Here denotes the quadrupole deformation of the nuclear wave function 
\P 7v 01 • The deformation of nuclear part becomes equal to /?o after the variation. It is noted 
that there are no constraint on the nuclear quadrupole deformation 7 and the deformation 
of the A single particle wave function. They have the optimum value after the variational 
calculation for each given value of (3 . 

Another constraint is on the A single particle wave function, 

V s = A\tp a )((p.\, (15) 
(r\cp 8 ) = exp{-pz/r 2 }, (16) 
v = yv x v y v z (17) 
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By using sufficiently large value for A, this potential forbids the A in s-wave. Therefore, by 
switching off and on this potential, we respectively obtain the hypernuclear state in which 
a A particle dominantly occupies in s- and p-waves. 
The total energy plus constraint potentials, 

E ' = hM^T + {v - } + ^ (18) 

is minimized using the frictional cooling method. The imaginary time development equations 
of the variational parameters are given as, 

dX n dE' 

It ~ hW*' [ ' 

X = Zi, z m ,ai, /3i,a m ,b m ,i' cn u i , (20) 

where \x is arbitrary negative real number. It is easy to proof that E' decreases as time 
develops, and after sufficient time steps we obtain the energy minimum under the constraint. 
By this method, we obtain the hypernuclear wave function for given f3 , the total parity and 
the A single particle orbital. In the present work, A dominantly occupies s- or p-wave and we 
shall use the notation A s and A p for them. Combined with the parity projection, we obtain 
four different configurations in which A s and A p couple to the positive- and negative-parity 
states of the core. They are denoted as (g) A s , £g> A s , (g) A p and ® A p in the 
following. 



III. RESULTS AND DISCUSSIONS 
A. general trend of the energy curves 

We have performed the variational calculation for ^Be, ^ 3 C, ^°Ne and ^Ne. To illustrate 
the change of nuclear deformation, Figure [1] shows energy curves of hypernuclear states 
with different configurations and corresponding normal nuclei as functions of deformation 
p. Each energy curve has an energy minimum shown by the open circle, and the binding 
energies, quadrupole deformation and radius at the minimum are listed in Table U The 
binding energy of A is defined as the difference of energy between the ground state of the 
core nucleus and the hypernuclear states, 

B A = B(i^X) - B{ A X g , s ) (21) 
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FIG. 1: (color online) Energy curves as function of nuclear quadrupole deformation /3 for (a) ^Be, 
(b)(c) ^ 3 C, (d)(e) ^°Ne and (f)(g) ^"Ne. (a),(b),(d) and (f) compare the positive-parity states of 
normal nuclei (^%) with the hypernuclear states of Vl/^-® A s and ^~^®A p configurations. (c),(e) and 
(f) compare the negative-parity states (^jy) with the hypernuclear states of VP^y® A s and \I'^ r ® A p 
configurations. Open circle shows the energy minimum on each curve. Energies of hypernuclei are 
shifted as shown in the figure for the sake of the presentation. 



12.0 



o > 10.0- 



&^ 2.0 

.3 
< 



0.0 



- — 


3 A s 


- 








- 


$ A p 





0.0 



quadrupole deformation p 



0.6 



FIG. 2: (color online) The single particle energy of A s and A p of ^ 3 C as function of the quadrupole 
deformation of core nucleus 12 C. Solid (dashed) line shows the energy of A s coupled to the positive 
(negative) parity state of 12 C. Dotted line shows the energy of A p coupled to the negative-parity 
state of 12 C. 

As the general trend, the shape of energy curve is not strongly modified by A particle 
except for ^ 3 C, and deformation (3 at the minima are slightly changed. In all cases, A s reduces 



quadrupole deformation. This is consistent with the cluster model calculations |22l . [27] 
and the (relativistic) mean- field calculations js-5] that demonstrated the reduction of /3 by 
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TABLE I: Calculated total and A binding energies B, B\ in MeV, quadrupole deformation pa- 
rameters j3 and 7, and the root mean square radius in fm at the minimum of each energy curve. 



Central values of observed energies 
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361 ] are also listed in parenthesis. The energies for 



A p states of ^Be and i C are estimated from the observed excitation energies given in Ref. 
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A s . On the other hand, it is found that A p increases 0. The magnitude of the change in 
quadrupole deformation is strongly dependent on the core nucleus, and the most prominent 
in ^ 3 C in which A s makes 12 C core spherical, while A p enhances the core deformation. The 
reason of the opposite behavior of A s and A p and the strong dependence on the core nucleus 
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is clearly seen in the single particle energy of A. Figure [2] shows the single particle energy of A 
(e\(f3)) in each parity and A single particle state in \ 3 C. Here e\((3) is defined as the difference 
between the binding energy of \ 3 C with the deformation and that of corresponding state 
of 12 C with the same deformation, 

e A (^) = Sis c (/3)-Si ac (/3). (22) 

It shows the Nilsson-model-like behavior of the A single particle energy. The binding of A s 
becomes shallower as deformation becomes larger. In the case of A p , A occupies the lowest 
p-wave that comes down as deformation becomes larger. Therefore, A s makes quadrupole 
deformation smaller and A p in the lowest p-wave makes it larger. The A single particle energy 
varies within a range of 1~2 MeV as function of quadrupole deformation that is smaller than 
the variation of the core nucleus energy. It is also case for other calculated hypernuclei. This 
explains why only A 3 C (FIG. CD (b)) manifests the drastic change in quadrupole deformation. 
Since the positive-parity state of 12 C is quite soft against the quadrupole deformation, small 
change in the A single particle energy can result in the large modification in quadrupole 
deformation. In other cases, change in A single particle energy cannot overcome much larger 
variation of the core energy and results in minor modification of quadrupole deformation. 
Therefore, we can conclude that the drastic change of quadrupole deformation by A particle 
occurs when the core nucleus is quite soft against quadrupole deformation within a range of 
2 MeV. Since the behavior of the energy curve is sensitive to the effective NN interaction 



371 ]. the drastic change in A 3 C may depend on the choice of NN interaction and it will 
be investigated in our future work. The behavior of the A single particle energy is also 
understood from the density distribution of the A particle and the core nucleus as shown in 
Figure [3J It shows that as nuclear deformation becomes larger the overlap between the A s 
and the core wave function becomes smaller (for example, compare £g> A s and Cg> A s 
of 20 Ne). It leads to the reduction of AN attraction. On the contrary, larger deformation 
makes the overlap larger in the case of A p and increases AN attraction (see Cg> A p and 
^® A p of A °Ne). In the case of A p , larger deformation reduces the kinetic energy that also 
contribute to the deeper binding of A p . 

Another issue to be mentioned is the reduction of the nuclear radius by A particle. In 
all cases, the radius of nuclear part is reduced, but the reduction (less than 5%) is much 
smaller than in the case of 7 A Li (20%) [lOj. More detailed discusstion will be made in our 
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next work. 



B. discussion on each hypernucleus 

The calculated total binding energies of 8 Be and ^Be underestimate the observed values by 
about 5MeV. The underestimation is common to all other hypernuclei. It will be resolved by 
performing the angular momentum projection and the generator coordinate method (GCM) 
that are usually performed in the study of normal nuclei by AMD. Indeed, in the case of 
20 Ne, it was shown that AMD calculation [15] reproduces the observed binding energy. The 
angular momentum projection and GCM will be performed in our next work. Despite of 
the underestimatin of the total binding energy B\ of A p is comparable with the observed 
value. 

The density distribution of A p in 9 Be <8> A p in FIG. [3]) clearly shows that this state 



corresponds to the supersymmetric (genuine) hypernuclear state [ljJ-UJjfl. The nuclear part 
has the pronouced 2a cluster structure and the A occupies the p orbital parallel to the 
symmetry axis. It is also interesting to note that A s state reduces the inter-cluster distance, 
while A p state increases it. 

^ 3 C manifests the drastic change in the quadrupole deformation. The A s makes 12 C core 
spherical, while A p state enhances deformation. It is noted that 12 C has the Op^/2 subshell 
closure configuration at small deformation and 3a cluster structure develops as deformation 
becomes larger. In other words, the nucleon spin is not satureted at small deformation, 



while that is alomost zero at larger deformation. The sophisticated AMD calculation [14 1 
has shown that the low-lying states of 12 C have mixed nature between the Ops/2 subshell 
closure configuration and 3a cluster structure, and the mixing strength is different for each 
state. Since A particle changes the deformation and spin property of 12 C, it will have 
influence on the AN spin-orbit splitting of \ 3 C 



38|-|40| 



Based on the cluster model calculation, the parity inversion in ^°Ne was suggested by 
Q 

Sakuda and Bando [7(. The core nucleus 19 Ne has the a+ 15 cluster state (J 7r =l/2j") 238 



keV above the ground state (J 7r =l/2 + ) that has (sd) 3 shell structure 41] . They concluded 
that A s coupled to the J 7r =l/2]~ state was more deeply bound than that coupled to the gound 
state, and the J n =l/2~ £g> A s configuration becomes the ground state of ^°Ne. They argued 
that the J n =l/2± state has dilute a+ 15 cluster structure and by the reduction of the 



10 



inter-cluster distance, A s gains larger binding energy than the J n =l/2^ <g> A s configuration. 
Our result shows the opposite trend to their result. Since the positive-parity state is more 
deformed than the negative-parity state, the binding of A s is weaker when it coupled to 
the negative-parity state. This trend is common to other calculations including the cluster 
model calculation for 21 Ne [22] . However, AMD fails to reproduce small excitation energy of 
the negative-parity state and it does not have prominent a+ 15 clustering, that are mainly 
due to the lack of the angular momentum projection and the GCM calculation. We will 
need more sophisticated AMD calculation to settle down this problem. 

The negative-parity state of 20 Ne has larger deformation than the positive-parity state. 
Therefore, the A s coupled to the positive-parity state is more deeply bound than that coupled 
to the negative-parity state. It is common to other hypernuclei studied here. On the 
contrary, A p is more deeply bound to the negative-parity state. Since number of nucleons 
in 20 Ne is large enough to bound A p , we can expect that the addition of A will generate 
a variety of bound rotational bands in fNe as disscussed by Yamada et al. 22)]. We will 
discuss ^Ne in detail in the forthcoming paper. 

IV. SUMMARY 

An extended version of AMD named HyperAMD has been introduced to investigate struc- 
ture of p-sd shell hypernuclei. The energy curves of ^Be, \ 3 C, /^Ne and ^Ne as functions of 
quadrupole deformation are studied. It has been found that A s reduces nuclear deformation, 
while A p increases it. It is due to the variation of the single particle energy of A as function 
of quadrupole deformation. The binding of A s decreases as deformation becomes larger, 
while that of A p increases. The variation of A single particle energy is within a range of 
1~2 MeV, that is rather small compared to the variation of the energy of the core nucleus. 
Therefore, the magnitude of the change of deformation strongly depends on the softness of 
the core nucleus against quadrupole deformation. Since 12 C is very soft against quadrupole 
deformation, it manfests the most prominent change of quadrupole deformation. This trend 
of the deformation change caused by A s and A p contradicts to the cluster model calculation 



for fNe 



but is consistent with other calculations. More sphisticated AMD calculation 



11 



will be needed to resolve this disagreement. 
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FIG. 3: (color online) Density plot of the intrinsic wave function at each minimum on the energy 
curve. Contour lines shows the density of the nuclear part and the color plot shows the A single 
particle orbital <p\. Dots in the figure shows the centroids of Gaussian wave packets on nuclear 
part Vl/jv- 
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